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Abstract. We describe quantum entropic functionals and outline a research program dealing 
with entropic fluctuations in non-equilibrium quantum statistical mechanics. 



1 Introduction 

Starting with papers [7, 8, 18, 19], the mathematical theory of non-equilibrium quantum statistical me- 
chanics has developed rapidly in the last decade. The initial developments concerned the theoiy of 
non-equilibrium steady states, the entropy production observable, and linear response theory (Green- 
Kubo formulas, Onsager reciprocity relations) for open systems driven by thermodynamical forces (say 
temperature differentials). This line of development was a natural direct quantization of the classical 
theory. In contrast, extensions of the classical fluctuation relations of Evans-Searles [h ] and Gallavotti- 
Cohen to the quantum domain have led to some surprises and novel classes of entropic functionals 
with somewhat striking mathematical stmcture and physical interpretation. 

A pedagogical introduction to our research program dealing with fluctuation theorems/relations in non- 
equilibrium quantum statistical mechanics can be found in [12]. This note can be viewed as a brief 
introduction to [ i ]. We sketch in telegraphic and simple terms the finite time/finite volume theory in the 
classical case (for comparison purposes) and in the quantum case, and comment on the resulting research 
program. The interested reader may consult [13, 12] for additional information. 
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2 Classical picture 

Consider a pair {E, (j)), where £ = {C,j}o<j<N is a finite phase space and cj) : £ ^ £ is a discrete time 
dynamics. For our purposes without loss of generality we may assume that 4>{Cj) = Cj+i (Ca^+i = Co)- 
Observables are functions / : <S — M and states are non-vanishing probability measures p on £. We write 
P{f) — Yli^ f{C)p{C)- Observables evolve in time as = / o 0*, t G Z, and states as pt{f) = p{ft)- 

The relative entropy of two states 

5(p,i.) = J]p(C)iogKC)MC)), 

c 

satisfies S{p, u) < and S{p, i^) = iff p = u. The Renyi relative entropy of order a € M is defined by 

5„(/,,z.) = iog^p(c)i-v(cr. 

C 

Our starting point is a dynamical system {£, 0, loq) where loq is a given reference state. We assume that 
uq is not constant (and hence not invariant under 0) and that wo(Ci) ~ ^o{CN-j)- This last assumption 
ensures that our dynamical system is time reversal invariant (TRI) with time reversal 0{(j) = (n-j- 

The entropy observable is 5o(C) = — logu;o(C)- The observable S* = (5^ — 5o)/t describes the mean 
entropy production rate over the time interval [0,t\. One easily verifies that for t > 0, 

L^o(S*) = -^5K,L^o)>0, (2.1) 

in accordance with the (finite time) second law of thermodynamics. The entropy production observable 
(or the phase space contraction rate) is defined by a{() = — log (a;i(C)/wo(C)) and satisfies S* = 

The relation (2.1) holds without the TRI assumption. The TRI however allows to refine the second law 
as follows. Let £tx = {( \ S*(C) = A} and p*(A) = ooo{£tx). An easy computation gives the celebrated 
Evans-Searles fluctuation relation 

p\-X) = e-V(A). (2.2) 

This relation implies a;o(S*) > and is saying more: the negative values of the mean entropy production 
rate ai^e exponentially suppressed in a universal manner. 

The classical entropic functional is defined by 

et(a) =logu;o(e-"*^*) . (2.3) 

The symmetry 

et{a) = et{l-a), (2.4) 

which holds for all a G M, is an equivalent formulation of the fluctuation relation (2.2). Clearly, 6^(0) = 
0, and hence et(l) = (this is sometimes called the Kawasaki identity [2]). The function a i-^ et{a) is 
convex, and ej(0) = — two(5^*)- The classical entropic functional satisfies 

et{a) = max S{p,ujq) — atp{Tj^), (2.5) 
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and 

et{a) = Sa{uJt,u;o). (2.6) 

The functional et{a) can be also described in terms of Ruelle transfer operators. For p € [1, oo[ we set 

ll/llp = Ecl/(C)P^o(C) and define 

Upit)f = f.te-^^" = /.^ei'^-'e"^'^". (2.7) 

Then 

Upih + 12) = Up{h)Up{t2), Up{-t)fUp{t)g = ftg, \\Up{t)f\\p = \\f\\p, 

i.e., Up is a group of isometrics of the space L^(£', wq) which implements the dynamics. In temis of this 
group, one has 

et(a) = log||f/p/,(t)l||^, (2.8) 

where 1(C) = 1- 

The results described in this section extend under minimal regularity assumptions to an essentially arbi- 
trary classical dynamical system [13]. 

3 Quantum picture 

Consider a pair (/C, H) where /C is a finite dimensional Hilbert space and H is a. Hamiltonian. Ob- 
servables are linear maps ^ : /C — )• /C (the identity map is denoted by 1) and states are strictly pos- 
itive density matrices p on /C. The observables evolve in time as At = e'*^Ae~^*^, t G M, and 
states as pt = e~'*^/9e'*^. We write p{A) = tr(Ap). The relative entropy of two states S{p, u) = 
tr(/3(logp — logz^)) satisfies S{p,i') < and S{p,v) = iff p = z^. The Renyi relative entropy is 
defined by 5q,(/9, v) = \ogii{p°'v^~"^). 

Our starting point is a quantum dynamical system (/C, if, wq), where is a given reference state. We 
assume that ujq does not commute with H and that the system is TRI, i.e., that there exists a complex 
conjugation on /C that commutes with H and wq- 

The entropy observable is Sq = — loga;o. S* = {St — So)/t is the mean entropy production rate 
observable and for t > the finite time second law holds: 

wo(S*) = -^5(^t,a;o) >0. (3.9) 

The entropy production observable (or the quantum phase space contraction rate) is o" = —i[H, logiLio] 
and S* = Jq CTsds. 

To aid the reader we describe one concrete physical setup. Consider two quantum dynamical systems 
{1Ciij.,Hii^,uiiI^), colloquially called the left and the right. Assume that initially the l/r system is in 
thermal equilibrium at inverse temperature Pij^, i.e., that ujij^ = e"'''/'"^'/'"/^^/^. Let IC = K.^® ICr, 
ujQ = uji (X) Ur, and H = Hi + + V , where V describes the interaction between the left and the right 
system. In this case 
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where = V] satisfies 

Jo 

and describes tiie energy flux out of tiie l/r system. 

We now turn to quantum fluctuation relations. TRI implies that the spectrum of S* is symmetric w.r.t. the 
origin. If Ptx is the spectral projection of S* onto A andp*(A) = ujo{Pxt), then the direct quantization of 
the fluctuation relation (2.2) is 

p\-X) = e-V(A). (3.10) 

Setting 

ei(a)=log.;o(e-"*^'), 

one easily verifies that (3.10) holds iff et{a) = et{l — a). However, one can show that the quantum 
Kawasaki identity et{l) = holds for all t iff [H,ujq] = and so the direct quantization of (2.2) fails. 
The standard route of observable quantization does not lead to quantum fluctuation relations. 

The following family of quantum entropic functional indexed by p G [1, oo] was introduced in [12]: 



epi(a) 



p 

1 — a o 2a c* 1 — d c — 

on — ot on 

g p "g p g p 



^ ' if 1 < p < oo, 



log tr 

log tr (g(l-a)'So+a-St) if p = OO. 

To motivate these functionals, note that 

Coot (a) = maxS'(;9,a;o) - atp(S*), 
p 

and so eoot(a) is the quantization of the variational formula (2.5). We also have 

e2t{a) = Sa{uJt,u}), 

and so 624 (a) is the quantization of (2.6). Regarding the other functionals, we need to introduce first the 
quantization of Ruelle's transfer operators. The quantization of the usual L^-norm is the Araki-Masuda 
L^-norm 

\\A\\P = tr(^\Aujl\^ 
The classical transfer operators (2.7) are quantized as 

Up{t)A = A.te^^''e~p^\ 

They satisfy 

Upih + 12) = Up{ti)Upit2), Up{-t)AUp{t)B = AtB, \\Upit)A\\p = \\A\\p, 
and for p G [1, oo[. 



ept{a) = log ||C/p/„(t)l| 
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Hence, the functionals ept{a) for p € [1, oo[ arise as the quantization of (2.8). 
The functionals ept{a) have the following properties. The symmetry 

Gptia) = ept{l - a), 

holds for all p and a and implies the quantum Kawasaki identity ept(l) = epj(O) = 0. The function 
a ep((a) is convex and the function p i— > ept{a) is continuous and decreasing. For all p, ep^(O) = 

-t6Jo(S*). 

The functional e2t{ct) has appeared previously in the literature in two seemingly unrelated contexts. 
In [I J] this functional was related to the fundamental concept of Full Counting Statistics (PCS) asso- 
ciated to the repeated measurement protocol of the entropy flow. Let Sq = \P\ be the spectral 
resolution of the entropy observable. The probability that a measurement of Sq at time t = (when 
the system is in the state wq) yields A is ujq{P\). After the measurement, the system is in the reduced 
state ujqP\/ujq{P\) which evolves in time as e~^^^ ujQP\e^^^ /ujq{P\). Denoting A' the outcome of a 
second measurement of Sq at a later time t > 0, the joint probability distribution of these two measure- 
ments is tr [e~^^^ P\e^^^ P\i) . It follows that the probability of observing a mean rate of entropy change 
= (A' - X)/t is 

X'-X=t<p 

The discrete probability measure is the PCS for the operationally defined entropy change over the 
time interval [0, t] as specified by the above measurement protocol. It is easy to verify that the e2t(a) is 
the cumulant generating function for the PCS, i.e., 

e2t(a) = log5^e-*"'^Pt(</.), 
<t> 

and the symmetry e2t(a) = e2t(l — a) yields the fluctuation relation 

In [20] the functional e2t(a) was motivated by the algebraic characterization of the Zubarev dynamical 
ensemble [21]. Let O be the vector space of all linear maps A : )C ^ JC equipped with the inner product 
{A, B) = iT{A*B). The relative modular- operator 

is a strictly positive operator on the Hilbert space (O, (•, •)). Let Qt be the spectral measure of -i log A^^j, 

1/2 

for the vector uj^ GO. Then 

e2t(a) = logj;e-"*'^gj(5), 

s 

and in particular, = Qt- This identification provides a striking link between the PCS and modular 
theory with far reaching implications. 
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4 Research program 

In telegraphic terms, our research program can be outlined as follows. 

(a) The first part of the program, carried out in [ iii], deals with development of the finite time theory 
of entropic functional in a general dynamical system setting. In the classical case this step is relatively 
easy [ ]. The quantum case is substantially more difficult and the full machinery of Tomita-Takesaki 
modular theory is required. The functional ept{a) and quantum Ruelle transfer operators are based on 
the Araki-Masuda theory of non-commutative L^-spaces [ ]. The step (a) can be viewed as an abstract 
unraveling of the mathematical structures underlying the entropic functional. These structures turn out 
to be of considerable conceptual and computational importance. 

(b) This step concerns the existence and regularity properties of the limiting functionals 

ep+{a) = lim \ept{a), (4.11) 

and is a difficult ergodic type problem that can be studied only in the context of concrete models. The 
existing results cover open spin-fermion and spin-boson systems (the proofs are based on the analysis of 
resonances of quantum Ruelle transfer operators, see also [ 1 6] for a pioneering work on the subject), and 
open locally interacting fermionic systems (the proofs are based on C*-scattering techniques). In some 
special cases (like the XY chain) the functionals ep^{a) can be expressed in closed form and analyzed 
in great detail. The step (b) of the program has been canied out in a series of papers and is a joint work 
with B. Landon, Y. Ogata, A. Panati, Y. Pautrat, and M. Westrich. It remains a challenge to extend these 
results to a wider class of models. We emphasize that the step (b) is meaningful only in the context of 
infinitely extended models (in other words, the thermodynamic limit must precede the large time limit). 
The thermodynamic limit is also needed for the physical interpretation of the finite time quantum en- 
tropic functionals of infinitely extended systems. 

(c) The Legendre transform of the limiting functional 62+ (a) is the rate function describing the large 
deviation fluctuations of the full counting statistics as t — > 00. In the linear regime and under suitable 
regularity assumptions the existence of 62+ (a) also implies the central limit theorem for the full count- 
ing statistics. For open quantum systems and in the linear- regime (neai^ equilibrium and for small a) all 
functionals ep+{a) reduce to Green-Kubo formulas for heat/charge currents. However, the development 
of linear response theory that goes beyond open quantum systems and allows for general thermody- 
namical/mechanical forces can be based only on the functional eoo+(«) (see [5, 13] for the results in 
the classical case that motivated this development). The functionals ep+(a), together with symmetries 
ep+(a) = ep+(l — a), can be viewed as extensions of the fluctuation-dissipation theorem to the far from 
equilibrium regime. 

(d) The functional 62+ (a) coincides with the Chemoff error exponent for the quantum hypothesis testing 
of the aiTow of time and this links quantum hypothesis testing, a rapidly developing branch of quantum 
information theory, to non-equilibrium statistical mechanics. This connection has been explored in [11]. 

(e) We have only discussed the entropic functionals defined with respect to the reference (initial) state 
of the system (they are sometimes called Evans-Searles type functionals). The Gallavotti-Cohen type 
entropic functionals are defined with respect to the non-equilibrium steady state (NESS) (the state that 
the infinitely extended system reaches in the large time limit). The Gallavotti-Cohen type function- 
als ai^e considerably more technical to introduce and study, and in the quantum case have a somewhat 
delicate physical interpretation. The Principle of Regular Entropic Fluctuations (PREF) introduced in 



6 



Entropic Functionals in Quantum Statistical Mechanics 



[13, 10] asserts that under normal conditions the limiting entropic functional of the Evans-Searles and 
Gallavoti-Cohen type are identical. Since under normal conditions the NESS is singular with respect to 
the reference state, the PREF can be viewed as a strong ergodic property of the physical model under 
consideration. 

(f) The developement of nonequilibrium statistical mechanics of open quantum systems in the Markovian 
approximation started with the pioneering work of Lebowitz and Spohn [15]. Since the Markovian ap- 
proximation is often the only technically accessible way to describe an open system, it is also important 
to develop a general theory of entropic fluctuations in this framework. The first attempts in this direction 
are due to Dereziriski, De Roeck and Maes [17, 3]. In [9], we derive fluctuation relations staiting from 
structural properties of the generator (Lindbladian) of the Markovian dynamics. 

5 Remarks 

Non-equilibrium statistical mechanics is a difficult subject and for many years our theoretical under- 
standing has been restricted to linear regime near equilibrium (linear response theory, the fluctuation- 
dissipation theorem). There ai^e good reasons for this: the richness and variety of non-equilibrium phe- 
nomena indicate that far from equilibrium physics may tolerate very few universal constraints. The 
recently discovered fluctuation relations of Evans-Searles and Gallavotti-Cohen are two universally valid 
constraints that hold far from equilibrium and reduce to linear response near equilibrium. These insights 
and subsequent developments (see [ ] for references) have dramatically altered our understanding of 
classical non-equilibrium statistical mechanics. The extensions of fluctuation relations to quantum do- 
main have led to further surprises that still remain to be fully explored. The research program outlined 
in this note is a first step in this direction. 
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